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Spin squeezing of an ensemble of N atoms in a high-Q cavity is shown to be enhanced by continuous 
measurement of photons in the cavity. A degree of spin squeezing proportional to iV -1 / 3 is attained 
as a quasi-steady state over a broad set of initial states of the atoms and photons. 
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I. INTRODUCTION 

Squeezed spin states (SSS) Q, in which quantum un- 
certainty of the collective spin in a direction orthogonal 
to the mean spin vector is suppressed below the standard 
quantum limit, have recently been realized experimen- 
tally 0jS0|- I n the SSS of atomic systems, entanglement 
between atoms is crucial for suppressing fluctuations in 
collective atomic spin, and therefore such systems are of 
potential interest in quantum information processing as 
well as in atomic interfcrometry. 

A number of methods to generate SSS in atomic sys- 
tems have been proposed, and these can be classified into 
several categories. A few examples that are relevant to 
the subject of the present paper are: (i) Squeezing is 
transferred between photons and atoms via their inter- 
actions 

H E H Using this property, SSS can be 
realized in an atomic system by copying the squeezing 
from the system of photons to the atoms. This scheme 
has been demonstrated experimentally 0. (ii) A probe 
laser beam propagating through atoms establishes en- 
tanglement between the atoms, with measurement of the 
probe laser induc ing sta te reduction in the atomic states, 
resulting in SSS 07113, E3l- 

This kind of measurement- 
induced-squeezing scheme has also been employed in ex- 
periments |3, 0|- (iii) The interaction between atoms and 
photons confined in a high-Q cavity entangles them, and 
states of both atoms and photons develop into squeezed 
states [1 0II1 E3. 

In the method (iii) proposed by Wineland et al. 
and by us H,S0, photons prepared in a coherent state 
and atoms prepared in a coherent spin state undergo in- 
teractions in a high-Q cavity, and after a certain interac- 
tion time SSS is realized temporarily, but not in a steady 
state. Precise control over the initial state and interac- 
tion time is therefore required. In the present paper, a 
novel scheme to generate SSS is proposed that combines 
methods (iii) and (ii) to circumvent the stringent require- 
ments of the precise control needed for method (iii) . SSS 
is thus obtained as a "quasi-steady state", regardless of 
the initial state, by performing continuous measurement 
of photons in the cavity. This drastically reduces experi- 
mental difficulties in producing SSS. Moreover, a degree 
of spin squeezing proportional to N^ 1 ^ , where N is the 
number of atoms, can be achieved using this method, in 



contrast to the N x / 4 attainable by method (iii) Q. 

This paper is organized as follows. SectionlTTldescribes 
the system we consider and develops a theory of continu- 
ous measurement of photon number in the cavity that is 
being pumped and contains atoms to be squeezed. Sec- 
tion HTT| shows that the atomic state in the cavity devel- 
ops into an SSS by continuous measurement of photons 
in the cavity. Section llVI discusses possible experimental 
schemes and concludes this paper. 



II. CONTINUOUS MEASUREMENT OF 
PHOTON NUMBER IN A PUMPED CAVITY 
WITH ATOMS 



A. Description of the system 



An operational theory of continuous photon counting 
has been proposed by Srinivas and Davies |l4j| . and mi- 
croscopic modeling of this theory has been reported in 
Ref. |l5| . The nonunitary state evolution due to mea- 
surement back action associated with the Srinivas-Davies 
model was investigated in Ref. [l^. In these studies, con- 
tinuous measurement of single-mode photons in a high-Q 
cavity was considered. We extend this model to the case 
of a pumped cavity containing N two-level atoms as il- 
lustrated in Fig. ^ 

We first consider the system without photon counting, 
namely, without the probe atoms and detector shown in 
Fig. n The high-Q cavity is assumed to sustain single- 
mode photons of frequency lu c and be pumped by a co- 
herent light source with an effective amplitude A and 
frequency uj v inside the cavity. The cavity contains two- 
level atoms with transition frequency w a , which couple to 
photons via the Jaynes-Cummings interaction |17| with 
coupling constant g. We assume that g is constant and 
the same for all the atoms, that is, we consider the Lamb- 
Dicke limit, where all atoms experience the same photon 
field. It is then convenient to introduce a set of collective 
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pumping 




cavity 



where <5 cp = lu c — aj p and 5 ap = aj a — tu p . 

We now consider continuous photon counting for the 
system described by the Hamiltonian Experimen- 
tally, photon counting may be performed using a pho- 
todetector placed in the immediate vicinity of the cavity, 
where photons leaking out of the cavity are detected. 
This kind of continuous photodetection can be simulated 
using a sequence of probe atoms that are injected through 
the cavity as illustrated in Fig.Q] The probe atoms 
are prepared in the ground state and injected into the 
cavity, with internal states detected upon leaving the cav- 
ity. If a probe atom is found to be in an excited state, 
we can say that one photon has been detected (one- count 
process). If it is found to have remained in the ground 
state, we can say that no photon has been detected (no- 
count process) |l5l Il6| . In the remaining part of this 
section we discuss nonunitary state evolution of the sys- 
tem in the one-count process, the no-count process, and 
more general processes. 



FIG. 1: A schematic of the proposed system. The TV two- level 
atoms interact with a single-mode photon field in a pumped 
cavity. The continuous measurement of photon number is 
simulated by a sequence of ground-state probe atoms injected 
through the cavity, and the state of each atom coming out of 
the cavity is measured by a detector. If a probe atom is found 
to be in an excited (ground) state, we assume that one photon 
(no photon) is detected. 



spin operators defined by 

JV 
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S, = ^E(|e J -)(e j |-|g i )(g. 

3=1 
1 N 



(1) 

(2) 



3=1 



where N is the number of atoms, and |gA and |ej) are 
the ground and excited states of the jth atom. We shall 
also use other collective spin operators defined by S x = 
(S+ + S-)/2 and S y = (S+ - S-)/(2i). 

The Hamiltonian for the system is given by 

H /h = u c h+uJ z +A(e l ^ t a+e- l ^ t a r )+g(aS++a*S-) 1 

where a (a T ) is the annihilation (creation) operator of the 
photon field, n = a* a is the photon-number operator, and 
A and g are taken to be real and positive without loss of 
generality. In Eq. ©, the first and the second terms on 
the right-hand side (rhs) describe the free Hamiltonian 
for the photon field and the collective atomic system, re- 
spectively; the third term describes the effect of coherent 
pumping of the photon field, and the last term describes 
the interaction between the photon field and the atomic 
system. In the rotating frame of reference e lu)p ( n+s ^\ip), 
the Hamiltonian @ is time- independent: 

H T ot /h = 5 cp h + 6 ap S z + A(a + tf)+g(aS + + tfS-), (4) 



B. State evolution in the one-count process 

If the passage time of each probe atom through the 
cavity is much shorter than the time scale of the state 
evolution described by Eq. the density operator of 
the entire system (i.e., cavity photons plus atoms exclud- 
ing the probe atoms) immediately after the one-count 
process at time t is given by (see Appendix [X] for the 
derivation) 



P(t+) 



ap(t)a' 



(5) 



where t+ denotes a time infinitesimally later than t and 
Tr the trace of the entire system. The denominator on the 
rhs of Eq. JSJ multiplied by 7 At gives the probability that 
a photon is detected between t and t + At, where 7 is the 
photodetection rate defined by Eq. I|A6() . The expression 
(J5J takes the same form as that of the Srinivas-Davies 
model 14]. Immediately after the one-count process (JSJ, 
the expectation value of an arbitrary operator O is given 
by 



Tr[p(t + )6] = (0) t+ = 



ah 



(6) 



When O is the photon-number operator h, we obtain 

((A«) 2 ) t 



(n)t+ = (n) t - 1 + 



(7) 



where An = n — (h). The one-count process thus de- 
creases the mean photon number by one and increases 
it by the Fano factor ((Ah) 2 ) t /(n) t of the photon field 
immediately before the one-count process |l6| . 

When O is some spin operator S, which commutes with 
a, Eq. © becomes 



(S)t + 



(n)t 



(S)t 



(AS At 



(8) 
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This indicates that upon the detection of a photon, the 
expectation value of the collective spin operator changes 
by an amount that depends on the correlation between 
the photons and atoms. As an example, let us consider 
an entangled state between an atom and a photon, where 
the ground and excited states of the atom are denoted 
by |g) and |e), and the vacuum and single-photon state 
by |0) and |1). Suppose that the state of the system be- 
fore measurement is given by (|g)|l) + |e)|0))/-\/2- Then 
(S z ) decreases from 1/2 to by the one-count process. 
This is because the measurement reveals that the cav- 
ity had contained a photon, and hence that the state of 
the system after the one-count process is given by |g)|0). 
Similarly, (S z ) increases upon the one-count process from 
1/2 to 1 for an entangled initial state (|g)|0) + |e)|l))/v / 2. 
In contrast, when the initial state is separable, i.e., 
(Sh)t = (S)t(n)t, the atomic quantities do not change 
by the one-count process. 



C. State evolution in the no-count process 

When no photon is detected between t and t + T, the 
density operator is shown to evolve according to the fol- 
lowing equation (see Appendix lAl for the derivation), 



P{t + T) = 



iH ac T/h 



P(t)e 



iHl T/h 



(9) 



where H nc = H r ot — ih/yh/2 is the non-Hermitian opera- 
tor that governs the state evolution during the no-count 
process. The denominator on the rhs of Eq. © gives 
the probability that no photon is detected between t and 
t + T. It follows from Eq. © that the master equation 
for the no-count process is given by 



am 

dt 



T\P(t)HL- H nc p(t)]+^(n) t p{t). 



(10) 



In fact, it can be shown by direct substitution of Eq. © 
into Eq. {TOJ that Eq. © is the solution of Eq. JE 



D. State evolution in the referring and 
non-referring measurement processes 

Time evolution under continuous photodetection is 
thus a stochastic process, in which the state evolution 
associated with one-count and no-count processes are de- 
scribed by Eqs. JSJ and 10, respectively. The probability 
of the one-count process occurring between t and t + dt is 
j{h)tdt and that of the no-count process is 1 — ~f{n)tdt. 
The number of photons detected and the times at which 
they are detected, therefore, differ run by run even for the 
same initial state. The measurement process in which 
we read the outcome of the measurement at every in- 
stant of time is referred to as referring measurement pro- 
cess ^(|. When the photons are detected at ti, ti, • ■ ■ ,t n 



from t = to t = T, the time evolution of the density 
operator is given by 

. (T) ^ P(T;h,--- ,t n )p(0)PHT;h,--- ,t n ) 



Tr[P(T; t u ■ ■ • , t n )p(0)Pt(T; t x , ■ ■ ■ , t n )} 



where 



P(T; ,t n ) = e 



iH ao {T-t n )/% SLe -iH ao {t n -t n - l )/% Si 

■ & e -iH™(t2-tl)/h &e -iHnetl/H^ ( 12 ) 



If the detector is switched on but we do not read out 
the measurement results, the state evolves into a statis- 
tical mixture of an ensemble of all possible outcomes, 
which is referred to as non-referring measurement pro- 
cess [lfj . The master equation describing the state evolu- 
tion in the non-referring measurement is obtained as fol- 
lows. In the non-referring measurement process, p(t + dt) 
is a statistical mixture of the density operator p±(t+dt) of 
the one-count process and the density operator po (t + dt) 
of the no-count process: 

p(t + dt) = ~f{h)tp\{t + dt)dt + (1 - j{h) t dt) p {t + dt). 

(13) 

Substituting pi(t + dt) for p(t + ) in Eq. (J5J gives 

-y(n) t p 1 (t + dt)dt = ~fap(t)a t dt. (14) 
Setting T = dt in Eq. 10 gives 



p (t+dt) = p(t)+-[p{t)Hl 



H nc p(t)}dt+l(n) t p(t)dt+0(dt 2 ). 

(15) 

Substituting Eqs. {H} and {TSJ into Eq. {T^Jl gives the 
master equation for the non-referring measurement: 



dp(t) 
dt 



h 



[2ap(t)a f - hp(t) - p(t)h] 



(16) 

When the cavity contains no atoms or when g = 0, 
we can show that the state of the photon field eventually 
evolves into a coherent state \a) with a = — A/(5 cp — 
i"f/2) in both referring and non-referring measurement 
processes, unless the initial state is orthogonal with \a*) 
(see Appendix lB|l . This fact is interesting in that the 
photon number never increases for the steady state \a) 
despite the fact that the cavity is being pumped and that 
no photons are detected. This can be understood from 
the equation describing the time evolution of the mean 
photon number in the no-count process 



d{ 



n) t 



dt 



= iA{ 



a + >t 



7 <(Ar 



(17) 



which can be obtained from Eq. (|10fl . For the steady 
state \a) with a = —A/(S cp — «7/2), the pumping term 
[the first term of the rhs of Eq. I|17fl ] balances with the 
last term in Eq. i|17|) which describes the decrease in the 
number of photons due to the state reduction by the no- 
count process. 
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III. GENERATION OF SQUEEZED SPIN STATE 
BY CONTINUOUS PHOTODETECTION 



A. Squeezed spin state 

Because the spin operators Sj (j — x,y,z) obey the 
commutation relation [S x , S y ] = iS x and the correspond- 
ing cyclic permutations, the spin vector fluctuates around 
the mean value (S) = ((S x ), (S y ), (S z )). For example, let 
us consider the lowest eigenstate of S z , |5, — 5), where 
5(5 + 1) is the eigenvalue of S 2 = 5 2 + S 2 + S 2 Z . The 
mean spin vector of this state is given by (S) = (0, 0, —5). 
Quantum fluctuations around this are isotropic and given 
by ((AS X ) 2 ) = ((ASy) 2 ) = 5/2, satisfying the minimum 
uncertainty relation ((AS x ) 2 )((AS y ) 2 ) = (S z ) 2 /4. If the 
fluctuation of one component perpendicular to the mean 
spin vector, say ((AS X ) 2 ), is suppressed below |(5 z )|/2, 
the spin state is defined as squeezed lj. 

A general definition of spin squeezing is obtained by 
rotating the above situation in spin space, when the mean 
spin vector (S) points in an arbitrary direction. The 
definition is thus given by Q] 



l(S)| 



(18) 



where S± is one spin component perpendicular to the 
mean spin vector, AS± = S±-{S±), and |(S)| 2 = (S x ) 2 + 
(S y ) 2 + (S z ) 2 . The squeezing parameter £ is unity in the 
coherent spin state 



n 



■e"" sm-| gj ; 



(19) 



The definition (|18|l of spin squeezing can also be applied 
to situations in which the spin state is a mixed state. 
The physical meaning of the SSS in atomic systems is 
discussed in Ref. 0. 

For an iV-atom system, there are N/2 + 1 (when N is 
even) or (N + l)/2 (when N is odd) subspaces that can 



be indexed by the value of the total spin 5 = N/2, N/2 — 
1, • • • ,0 (1/2). In the present paper, we restrict ourselves 
to the subspace of 5 = N/2. This corresponds to the case 
in which the initial state is the coherent spin state, which 
is easily prepared in experiments. Since the Hamiltonian 
ifljl commutes with S 2 , the subspaces do not mix by time 
evolution. 



B. Analytic results 

For simplicity, we consider the case of S cp = 6 ap = 0, 
i.e., the case in which both the coherent pumping field 
and the atomic transition energy are resonant with the 
cavity mode. The non-Hermitian operator H nc describ- 
ing the no-count process is then given by 



H nc /H = A(a + a r ) + g(aS+ + a f 5_) - yri 



(20) 



where we define a = g5_ + A. 

As will be shown in Sec. IIII CI the generation of SSS 
is effective for j/g > 5, and hence we derive an approx- 
imate master equation of the no-count process for large 
7. In the following derivation, we assume 7 > yl and 



A - gS. Defining p'{t) = e^ th / 2 p(t)e~< m / 2 , 
the master equation 11(J|) for the no-count process as 



we rewrite 



dp'jt) 
dt 



ip 



~yt . -I. . 

e 2 a'a 



(e~^aa^ + e^a)a*j p'(t) + j{h) t p'(t). 

(21) 



Since 7 > 4, the photon field is almost always in the 
vacuum state, which is justified below. We may there- 
fore approximate p'(t) as p'(t) ~ |0)(0| <£> p s (t), where 
p s (t) denotes the reduced density operator for the spin. 
Integrating Eq. (|21|1 iteratively gives 



p(t + At) = 



r rt+At 

<g> <^ p s (t) +7 / dti(n) tl p s {t) 

dt 2 dt x e ^ 2 (a^dp s {t)+ p s {t)d^&) + 1 2 (h) tl (h) t2 
t Jt L 

1)(1|<E> / dt 2 dt^e 1 ^ 2 ap s (t)a f + ••• , 



(22) 



where we have kept only relevant terms. We take At such that 'yAt ^> 1 and 7AtA 2 /7 2 -C 1 are satisfied. From the 
last term of Eq. (|22|l . we find that the mean photon number is given by (n)t+At — ^1~ 2 &)t, which is 0(A 2 /j 2 ). 
The last term and the term including r Y 2 (n)t 1 (fi}t 2 m Eq. I|22|) are therefore 0(A 2 /j 2 ) and can be ignored. Thus, 
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Eq. 1221) can be approximated as 
P(t + At) 



<^ p 8 (t) [0-Vp s (t) +p s (t)<T%] + (CTV) t p s (t) 

I 7 7 



(23) 



Since, the last two terms in the curly brackets in Eq. 1(23(1 
are 0(7 At A 2 /j 2 ), which is assumed to be small, we ob- 
tain an effective master equation of the spin state for the 
no-count process: 

= -^&p a (t)+p&)tfa\ + -{tfa) t p&). (24) 
ot 7 7 

It follows from Eq. 1(24(1 that the system evolves into 
the lowest eigenstate of the operator 



o'o = g 2 S+S- + gA(S+ + S- 



A 



(25) 



if the initial spin state is not orthogonal with the lowest 
eigenstate. To see this, let us expand the density operator 
as 



P{t) = c„ m |</>„) (0„ 



(26) 



where \4> n ) is the eigenstate of the operator ((25(1 with 
eigenvalue e n , i.e., a^a\4> n ) = e n \cf> n ). The master equa- 
tion i(24|) is rewritten as 



c nm (t) = -(2(o- t a) t - £„ - e rn )c nm (t), (27) 
7 



which can be solved to give 



c n m(t) = exp 



d£ 



7 



. . (28) 

Thus, the system develops into the state that minimizes 
£n+£m, that is, the lowest eigenstate of the operator 1)25(1 . 
if the corresponding initial coefficient c nm (0) is not zero. 
We note that Eq. (|2"5)l is independent of 7, and therefore 
small experimental fluctuations in 7 do not affect the 
steady state. Since the relaxation coefficient in Eq. I|24() 
is proportional to 7 -1 , the relaxation becomes slower for 
larger 7. This is because the number of photons in the 
cavity is smaller for larger 7, and therefore the state re- 
duction caused by the no-count process becomes less ef- 
fective. 

The steady state of the master equation ((2411 is the SSS, 
as will be numerically shown in Sec. IIII CI We evaluate 
perturbatively the degree of spin squeezing (|TH|l of the 
steady state. When i > j, the first-order perturbation 
with respect to g/A yields (Appendix Q 



C = i-^(iv-i) 



< 1, 



(29) 



showing that the state is an SSS. When A <C g, the 
second-order perturbation with respect to A/g yields 



(Appendix 

2A 2 

showing that the state is again an SSS. 

C. Exact numerical results 

We numerically simulate the stochastic process of the 
continuous measurement described by Eqs. (0 and i|l(J|) 
with H nc given by Eq. (|20|) . We divide the time into 
a small interval At such that r y(n)tAt is always much 
smaller than unity. We then generate a random number 
r distributed uniformly between and 1 for each time 
interval. If < r < 7(n)tAi, we assume that the one- 
count process occurs between t and t + At and calculate 
the density operator at t + At according to Eq. ((SJ. If 
7(n)tAi < r < 1, we assume that the no-count process 
occurs during t and t + At and calculate the density op- 
erator at t + At according to Eq. lfTU|) . 

Figure|21shows an example of the stochastic time evolu- 
tions of various quantities in the referring measurement 
process with A/g = 5, j/g = 10, and N = 10, where 
the photon field is initially in the vacuum state and the 
atoms are in \S, — S). The impulses that can be seen at 
the bottom of Fig. [3(c) indicate the times at which pho- 
tons are detected. They are clearly bunched, reflecting 
the fact that g^ = (a) 2 a 2 )/(h) 2 in the cavity is larger 
than unity most of the time as shown in Fig. [21 (a). The 
bunching nature of photons makes the duration of each 
no-count process longer, helping the system reach the no- 
count steady state (shown as plateaus in Fig.EJ. We note 
that in the no-count steady state, the SSS is generated 
with ( ~ 0.56 [solid curve in Fig. [21(c)], where ( is defined 
in Eq. JTSJ. In the non-referring measurement process, 
in contrast, no SSS is obtained [dashed line in Fig. 0(c), 
obtained by numerically integrating Eq. p6|l ]. indicat- 
ing that the information concerning the detection times 
of photons is crucial for obtaining the SSS. It should be 
emphasized that in such a referring measurement process 
£ is smaller than unity most of the time, and therefore 
the SSS is obtained as a quasi-steady state. 

Figure [3| (a) shows the quasi-probability distribution 
(QPD) corresponding to the spin state at t ~ 15, when 
the no-count steady state is reached. We can see that 
the spin fluctuation is squeezed in the azimuthal direc- 
tion. However, at t = 16 the spin state is disturbed by 
the photon counts and the corresponding QPD shown in 
Fig. [31 (b) exhibits a pattern similar to a rotated Dicke 
state. 
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FIG. 2: The solid curves show examples of time evolutions 
of (n), g< 2 ' = (a t2 a 2 )/<n), {§,}, [<SS> ], and C as defined in 
Eq. 1181 under continuous photodetection with N = 10, A = 
5g, and 7 = lOg. The initial state of the photon field is 
the vacuum state, and that of the atoms is \S,S Z = —S). 
The impulses in (c) indicate the times at which photons are 
detected. The dashed lines show the values for the steady 
state in the non-referring measurement. 



(a) t= 15 




FIG. 3: Gray-scale images of the quasi-probability distri- 
butions of the spin state (9,cj>\p 3 \6,<f>}, where p a is the re- 
duced density operator of spin and \9, <j)) is the coherent spin 
state defined in Eq. 1191 . (a) At t = 15 in Fig. [5] (no-count 
steady state) and (b) at t — 16 (just after a few photons were 
counted) . 



TABLE I: Power laws of the optimal squeezing parameter ( 
with respect to the number of atoms N for various squeezing 
schemes. 



c 



two-axis twisting fl] 
one-axis twisting fl] 
Jaynes-Cummings (with measurement) 
Jaynes-Cummings (without measurement) 



iV" 1 

JV-2/3 
N -l/3 
JV-1/4 



rizes the dependence of C on N in these schemes. 



The degree of spin squeezing, £, depends on the pa- 
rameters A, 7, and N. Figure 0] shows the dependence 
of Cst on A and 7 for N — 2 and 10, where Cst is the 
steady-state value of C in the no-count process. We see 
that the SSS is obtained for a wide range of parameters. 
We find that in both cases of TV = 2 and TV = 10 the max- 
imum squeezing is obtained for A/g ~ N/2 = S, and we 
have numerically confirmed that this finding holds true 
for other values of N. For 7 ^> A, Cst depends little on 7, 
in agreement with the analytic results of Eqs. I|25|) . 1)29(1 . 
and PU||, 

Figure shows the dependence of Cst and the angle of 
the spin vector 9 s t on N, where the squares are obtained 
by solving the master equation (|10f) and the circles by 
diagonalizing Eq. (|25|l . Both plots are in excellent agree- 
ment, confirming the validity of the analytical treatment 
in Sec. IIII Bl It is interesting to note that both Cst and 
6 s t are proportional to TV -1 / 3 . The other models of spin 
squeezing also exhibit such power laws: TV" 2 / 3 for one- 
axis twisting and TV -1 for two-axis counter twisting 0. 
In a high-Q cavity, the Jaynes-Cummings interaction be- 
tween photons and atoms generates spin squeezing pro- 
portional to TV" 1 / 4 0, demonstrating that the present 
scheme is more effective for larger N. Table I summa- 



IV. DISCUSSIONS AND CONCLUSIONS 

When the number of excited atoms is much smaller 
than the total number of atoms N, the Holstein- 
Primakoff transformation [Tsj can be approximated by 
S- ~ {2S) l ' 2 b and S+ ~ (2S) 1 / 2 b\ where b and fct 
are the bose operators satisfying [b,w] = 1. The non- 
Hermitian operator describing the no-count process H2U|) 
then becomes 

H nc /h = A(a + a 1 ) + G(aP + a r b) - yn, (31) 

where G = {2S) 1 ^ 2 g. In a similar manner to the deriva- 
tion given in Appendix El one can show that during 
the no-count process the system develops into |0) a |/3)b, 
where |0) o is the vacuum state of the photon field and 
\f3)b = e 13 ^ 1 0) b is the coherent state of the boson with 
f3 = —A/G. Therefore, the bosonization of spin elim- 
inates the squeezing properties of atoms, and thus the 
spin algebra, which plays a role in nonlinear saturation, 
is crucial for obtaining squeezing in the present scheme. 

Experimentally, the model presented may be imple- 
mented in various systems, such as a microwave cavity 
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FIG. 4: The dependence of the squeezing parameter £ s t on A 
and 7 for the no-count steady state with (a) TV = 2 and (b) 
TV = 10. The legends on the right show the gray scale of the 
value of £st. 



with Rydberg atoms [l!j, an optical cavity with cold 
atoms |2Cj . and trapped ions coupled with their center- 
of-mass motion [T^ . For example, let us consider a sys- 
tem in which Rb atoms in the circular Rydberg state 
pass through a superconducting microwave cavity [l9| . 
The coupling constant between the atoms and the cav- 
ity photons can be varied from to 24 kHz [l^. If we 
take k/(2tt) = 2 kHz, g/(2n) = 100 Hz, r p = 10 /its, and 



TV = 10, the conditions (|A3I) are satisfied clS /^^"p ~ 0.1 
and gSYp — Ar p ~ 0.03, where k is the coupling con- 
stant between probe atoms and photons, and t p is the 
passage time of a probe atom through the cavity (see Ap- 
pendix [SJ. These parameters give j/g = k 2 t„/ q ~ 2.5, 
and the lifetime of the Rydberg state ~ 30 ms [lj| im- 
poses the constraint gt < 19. The cavity relaxation time 
now reaches 0.3 s |2l|. Thus the scheme is experimentally 
feasible. 

In conclusion, we have studied the system in which 
TV two-level atoms interact with photons in a pumped 
cavity, and have shown that the atomic state develops 
into an SSS by continuous measurement of photon num- 
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FIG. 5: The TV dependence of the minimum squeez- 
ing parameter f s t and the corresponding spin angle 8 st = 
tan^ 1 ({S z ) I (Sx)) (inset) of the no-count steady state, where 
(S x ) < 0, (S y ) = 0, and (S z ) < 0. Squares were obtained by 
solving the master equation I I Oil and circles by diagonalizing 
Eq. (25) . 



ber in the cavity. In particular, the SSS is obtained as 
a quasi-steady state during the no-count process. The 
present method does not require fine tuning of the ini- 
tial state and the time and strength of interaction, and 
therefore significantly facilitates experimental realization 
of the SSS. 
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APPENDIX A: A MODEL OF THE 
CONTINUOUS MEASUREMENT OF PHOTON 
NUMBER 

For continuous measurement of photon number, the 
Srinivas-Davies "ansatz" |14| is usually employed. How- 
ever, it is by no means obvious whether it is valid in 
our system, and therefore, we explicitly derive nonuni- 
tary time evolutions JS) and JH) using the continuous 
photodetection model of Ref . [l5j ■ 

We consider a situation in which probe atoms prepared 
in the ground state |g) pass through the cavity succes- 
sively, and the internal states of the atoms coming out 
of the cavity are detected by, e.g., an ionization detec- 
tor (Fig.nj. When a probe atom is detected to be in the 
excited state |e), we regard this event as a one-count pro- 
cess. When the atom coming out of the cavity remains 
in the ground state, we regard this event as a no-count 
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process. 

The time evolution as the probe atom passes through 
the cavity is described by the Hamiltonian 

H p = H* ot + hn(a& + + tf&J), (Al) 

where H™ 1 is given in Eq. Q , <r + and <7_ are the raising 
and lowering operators of the probe atom, and k is a cou- 
pling constant which we assume to be real. We assume 
that the transition frequency of the probe atom is reso- 
nant with the cavity mode. When the probe atom enters 
the cavity at time t and comes out of it at t + t p , the 
density operator p sp for the system and the probe atom 
evolves as 



When successive probe atoms are found to be in the 
ground state over a time interval T, the density oper- 
ator of the system after such a no-count process is given 
by Eq. ©. 

Probe atoms were introduced to explicitly derive 
Eqs. JSJ and 10. However, these expressions may also 
describe a more general system under continuous pho- 
todetection, such as the situation in which photons leak- 
ing out of the cavity are detected by a photodetector. 



APPENDIX B: STEADY STATE IN A PUMPED 
CAVITY UNDER CONTINUOUS 
PHOTODETECTION 



p sp (t + r p ) = flt)®|g><g|--^[#p,/3(t)®|g)<g|] 



2h 



M#p[#p,p(t)®|g>(i 



(A2) 



where p(t) is the density operator for the system (exclud- 
ing the probe atom). 

Here we assume that r p is sufficiently small so that 
third and higher-order terms in Eq. (|A2() can be ignored. 
This approximation is valid if the change in the state of 
the system is small during r p . This condition is satisfied 
if the following conditions are met: 

At' 1 < 1, gSr' 1 < 1, S^t' 1 < 1, ^r" 1 < 1. 

(A3) 

When the probe atom is detected to be in the excited 
state |e), the density operator of the system is projected 
into 

Tr pro bc[|e)(e| j o sp (i + r p )] = K 2 T*ap(t)a f . (A4) 



The normalized form of Eq. l|A4jl gives the density oper- 
ator immediately after the no-count process: 



ap(t)a) 
Tr[ap(t)tf] ' 



(A5) 



where t + denotes a time infmitesimally later than t. The 
probability that the one-count process l|A4(l occurs is 
given by the trace of Eq. (|A4(1 . i.e., 7T p (n) t , where 



7 = K t p 



(A6) 



can be interpreted as the rate of photodetection per one 
photon. 

When the probe atom is detected in its ground state 
\g), the density operator of the system is projected into 

Tr probc [|g) (g\Pep(t + r p )] 

9 

IT A T „ 

~ p{ t) -f[HS*,m] - ^ H o 0t > [H T o ot ,P(t)] 



[np(t) + p(t)n] 



(A7) 



We show that the state of the system described by the 
Hamiltonian (0} eventually evolves into a coherent state 
in both referring and non-referring measurements. 

The non-Hcrmitian operator that describes the no- 
count process (0 is rewritten in the rotating frame by 

H nc /h = A(a + a f ) + (S cp - ij/2)h 

-«( 4,+ iOK)-f 

= iiefo^neSW)-^, (Bl) 

where we define f2 = S cp — i-f/2. We then obtain 

e -m no t/n K e 4(a+at) e -intft e -^(a+at) 

OO 

\^ p -m5o P t-n7*/2 e i|(a+at)| n w n | e -4(a+a t )_ 



n=0 



(B2) 



For jt ^> 1, the term corresponding to n = 0, i.e., 



^(a+a+) 



■^(a+a f ) | a ^Q,*| becomes dominant in 
the summation, where \a) is the coherent state with 
a = — A/Vi. We thus obtain 



-iH na t/h\ 



lim ■ 

t — >oo 



(B3) 



for (a*\ip) 7^ 0. Therefore, during the no-count process 
(JSJ the state develops into \a) unless the initial state is 
orthogonal to \a*). If (a*\ip) — 0, the term of n = 1 in 
the summation of Eq. ljB2fl becomes dominant, and the 
state develops into a displaced number state. 

If the photons are detected at times ti,t%, - • • , t n , the 
time development operator (|12|) for the referring mea- 
surement takes the form 



P(T; ti, • • ■ , t n ) 

e -iA M T/» nJ j =i 



— iQtk - 



.4 



— iQtk 



From Eqs. (|B3|) and (|B4|) we obtain 

lim P(T;ti, ■ ■ ■ ,t n )\ip) cx |a). 



.(B4) 



(B5) 
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Thus, the state develops to |a) in the referring measure- with the squeezing parameter calculated to be 
ment. This result is also true for non-referring measure- 
ment, since Eq. l|B5jl depends neither on the number of 

detected photons nor on the times at which photons are C = 1 — -^—(2S — 1) (C3) 

detected. AA y >' { ' 



APPENDIX C: PROOF OF EQS. AND (l30ll 

We calculate perturbatively (2^ the squeezing param- 
eter for the lowest eigenstate of the operator 1|25[1 . 

When g/A <C 1 , the first term on the rhs of Eq. l|2*5|l 
can be treated as a perturbation. For convenience, we 
rotate in the spin space by 7r/2 with respect to the S y 
axis. The operator i|25|) then becomes 



1 



2 A 

A 



S z 



A 2 



(S 2 



S x 



(CI) 



The unperturbed ground state is then taken to be the 
eigenstate of S z with eigenvalue — S, which we write as 
\S, — S) . Performing first-order perturbation with respect 
to the last term of Eq. IjClf) . we obtain the ground state 
as 



\S,-S) + 



2Sg 



AA 



1 



S, 1 - S) + -y/2(2S - 1)\S, 2~s)y 

(C2) 



When A -C g, the second term on the rhs of Eq. lt2*5|) 

may be treated as a perturbation. In this case, we must 
perform second-order perturbation to obtain spin squeez- 
ing. The ground state is then given by 



1 - 



.4 2 



ASg 2 
A 2 1 



\S,-S) 



g 2 v/45(25- 1) 



2Sg 
S,2-S), 



\S, 1 - S) 



(C4) 



with the squeezing parameter calculated to be 



C = i 



.4 2 



25V' 



(C5) 
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